function [dzt,zl,it] = cpleja(it,ntotal,nleja,zt,zl,dzt)
% From the it previous Leja points calculate the it+1 point.
      
   
      dmax=0.;
      for jj=1:ntotal

% Distances.
         if (it==1) 
            dzt(jj)=abs(zt(jj)-zl(it));
         else
            dzt(jj)=dzt(jj)*abs(zt(jj)-zl(it));
         end

% Maximum.
         if (dzt(jj)>dmax) 
            dmax=dzt(jj);
            kr=jj;   
         end 
      end 

% Next Leja point.
      it=it+1;
      zl(it)=zt(kr);

      end
